Weakly bound atomic trimers in ultracold traps 
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The experimental three-atom recombination coefficients of the atomic states 23 Na|_F = 1, mi? = 
— 1), 8r Rb|F = = —1} and 85 Rb|F = 2, iuf = —2), together with the corresponding two-body 

scattering lengths, allow predictions of the trimer bound state energies for such systems in a trap. 
The recombination parameter is given as a function of the weakly bound trimer energies, which are 
in the interval 1 < m(a/h) 2 E3 < 6.9 for large positive scattering lengths, a. The contribution of a 
deep-bound state to our prediction, in the case of 85 Rb|_F = 2,mp = —2), for a particular trap, is 
shown to be relatively small. 



The formation of molecules in ultracold atomic traps 
offers new and exciting possibilities to study the dy- 
namics of condensates [1]. It was reported the forma- 
tion of Rubidium molecules 87 Rb2 in a bosonic con- 
densate, which allowed to measure its binding energy 
with unprecedented accuracy [2]. Ultracold Sodium 
molecules 23 Na2 have also been formed through photoas- 
sociation [3]. However, nothing has been reported till 
now about formation of molecular trimers in cold traps. 
The first information one is led to ask is the magnitude of 
the binding energy of trimers in a cold trap. Two-body 
scattering lengths of trapped atoms are well known in 
several cases, as well as their closely related dimer bind- 
ing energy. In the limit of large scattering lengths, it 
is necessary to know in addition, one low-energy three- 
body observable to predict any other one. In this case, 
the detailed form of the two-body interaction is not im- 
portant [4,5]. The recombination rate of three atoms in 
the ultracold limit, measured by atomic losses in trapped 
condensed systems, can supply the necessary informa- 
tion to estimate the trimer binding energy. For short 
range interactions, the magnitude of the recombination 
rate of three atoms is mainly determined by the two- 
body scattering length, a [6]. However, it is important 
to remark that, still remains a dependence on one typical 
low-energy three-body scale [4,5]. Indeed, it is gratifying 
to note that all the works on three-body recombination, 
consistently, present a dependence on a three-body pa- 
rameter in addition to the scattering length [7-9]. 

The aim of the present work is to report on how one 
can obtain the trimer binding energy of a trapped atomic 
system, from the three-body recombination rate and the 
corresponding two-body scattering length. For this pur- 
pose, we use a scale independent approach valid in the 
limit of large positive scattering lengths (or when the 
interaction range goes to zero), obtained from a renor- 
malized zero-range three-body theory [4], which relates 



the recombination rate, the scattering length and the 
trimer binding energy. Considering the experimental val- 
ues of the recombination rates and scattering lengths 
given in Refs. [10-13], the method is applied to predict 
the trimer binding energies of 23 Na|_F = l,mj? = — 1), 
87 Rb|.F = l,m F = -1), and 85 Rb|F = 2,m F = -2), 
where \F,mp) is the respective hyperfine states of the 
total spin F. We note that the bound-states considered 
here are in fact high- lying resonances, not true bound 
states, as they can decay into lower-lying channels. 

The validity of our approach is restricted to sufficiently 
diluted gases, because all the scaling relations are de- 
rived for three isolated particles. Also, when the scat- 
tering length is tunned via external field in a trap, the 
parameters are different from the vacuum values, and 
consequently our predictions only apply to that particu- 
lar experimental conditions. For the trapped gases that 
we are analyzing, the diluteness parameter pa 3 (where p 
is the gas density) should not be much larger than one, 
otherwise one needs to consider higher order correlations 
between the particles. Indeed, we observe that, in gen- 
eral, for the analyzed condensed systems, the diluteness 
parameter is much smaller than one. Even in the case of 
85 Rb, where the considered scattering length is obtained 
via Feshbach resonance techniques [13], the diluteness 
parameter is about 1/2. 

Another relevant remark, pointed out in Ref. [14], is 
that the recombination into deep bound states can affect 
the theoretical results that are based on calculation of 
this rate into shallow states alone. This additive contri- 
bution depends on one more constant, beyond the three- 
body scale. However, the fitted contribution of the re- 
combination into deep bound states, is fortunately much 
smaller than the contribution of the shallow bound state, 
as found in the case of 85 Rb [14]. Such evidence supports 
our estimatives of trimer energies, when a > 0, that are 
obtained by only considering the contribution of recom- 
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bination rates into the shallow state. 

The values of a are usually defined as large in respect 
to the effective range r , such that a/r >> 1. The low- 
energy three-boson system presents, in this limit, the 
Efimov effect [15], where an infinite number of weakly 
bound three-body states appears. The size of such states 
are much larger than the effective range. The limit 
a/r — > oo can be realized either by a — > oo with ro 
kept constant or by ro — > with a constant. In the last 
case, the limit of a zero-range interaction, corresponds to 
the Thomas bound-state collapse [16]. In this respect, 
the Efimov and Thomas limits are equivalent; or, dif- 
ferent aspects of the same physics [17]. The Thomas- 
Efimov connection is also reviewed in Ref. [5]. In the 
limit a/ro — > oo, the details of the interaction for the 
low-energy three-body system are contained in one typi- 
cal three-body scale and the two-body scattering length 
(or the dimer bound-state energy, E 2 ); they are enough 
to determine all three-body observables [4]. Consider- 
ing, for example, the trimer binding energy (E 3 ) as the 
three-body scale, any three-body observable (0 3 ) that 
has dimension of [energy]' 3 , in the limit of r — > 0, can 
be expressed as 



O3 = E^ 2 (E 2 /E 3 ) = E^ 3 (E 2 /E 3 ) . 



(1) 



The dimensional factor in front of the above equation 
(1) is chosen for convenience as E 2 or E 3 . The scaling 
function in each case is T 2 or T 3 . The existence of the 
scaling limit for zero range interactions was verified in 
Refs. [18,19]. In practice, such limit is approached by 
the excited state of the atomic trimer obtained in re- 
alistic calculations, allowing as well the theoretical in- 
terpretation of those excited states as Efimov states [19]. 
Here, we observe that the binding energy E 3 refers to the 
magnitude of the total energy of the bound-system; the 
binding energy with respect to the two-body threshold is 
defined as S 3 = E 3 - E 2 . 

The rate of three free bosons to recombine, forming a 
dimer and one remaining particle, is given in the limit of 
zero energy, by the recombination coefficient [5,6] 
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K 3 = — a a, 
m 



where a is a dimensionless parameter and m the mass of 
the atom. When a > 0, the recombination parameter a 
oscillates between zero and a maximum value, which is 
a function of a, as shown in Refs. [7] ( a < 68.4 ), [8] 
(a < 65) and [9] (a < 67.9). With amplitude a max and 
phase S, we can write it as [5] 



(1.01 ln(a) + £) 



(3) 



where 6 depends on the interaction at short distances. 
The physics at short distances, in the three-boson sys- 
tem, is parametrized by one typical three-body scale, 
which we have chosen as the unknown trimer binding 
energy. So, by using the general scaling given by Eq. (1), 



one can explicitly express the functional dependence of 
aasoEa (^/ E 2 / E 3 ^j , considering that, for large scat- 
tering lengths we have 1/a — y mE 2 /h 2 . To exemplify 
the scaling form of a, we can rewrite Eq.(3) such that 
the \J E 2 jE 3 dependence is explicit. Therefore, 



a = a m axSin — 1.01 In i 



+ A 




(4) 



where A (y£ 2 /£ 3 ) = S - 1.01 In (^JmE 3 /h 2 ^j . Our 

next task is the calculation of the scaling function, by 
using the renormalized subtracted Faddeev equations [4] . 

The three-boson recombination coefficient at zero- 
energy is derived from the Fermi's golden-rule as 
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where 7^f is the transition amplitude between the initial 
(i) and final (f) momentum states, which are normalized 
as plane-waves: (r\p) = exp(— i(p//i).r) /(2ttH) 3 ^ 2 . The 
number of atoms N in the condensed state decreases, due 
to the recombination process, as: 



1 dN 3 2 



(6) 



For each recombination process three atoms are lost, jus- 
tifying the factor 3 in the numerator. The factor 3! in the 
denominator appears only in case of condensed systems; 
it counts for the number of triples in such state [11]. 

Considering the symmetrized scattering wave- function 
for the initial state of three free particles, |$o) = 
(1/V3) J2i=i \QiiPi)y we obtain the transition amplitude, 
in terms of the Faddeev components of the three body 
T-matrix, T 4 (£), as 



(k l ^ k) \[T J (E) + T k (Em ) 



(7) 



where (i,j,k) = (1,2,3) and cyclic permutations, and E 
is the energy of the scattering state is E = 3qf / (4m) + 
(2) Pi/m = 3A: 2 /(4m) — E 2 . <fj is the Jacobi relative mo- 
mentum of the particle i in respect to the center of mass 
of particles j and k, Pi is the relative momentum of j 
in respect to k, and fcj is the relative momentum of the 



free particle i in the final state. \^ k ^} is the normalized 
two-body bound state wave function of the pair (jk). The 
calculation of the Faddeev components is performed with 
the use of the subtracted approach given in Ref. [4] , such 
that 

Ti(E) = U (e - {1 + ( G +(E) G (E,)) x 

x(T (E)+T k (E))} , (8) 

where E^ = —/j, 2 /m is the subtraction energy scale with 
/i a constant in momentum units. It is possible to vary 
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ji without changing the physics of the theory as long as 
the inhomogeneous term of Eq. (8) is modified accord- 
ing to the renormalization group equations [20] . U is the 
two-body t— matrix for the subsystem of particles (jk). 
For E = and zero-range potential the corresponding 
matrix elements are given by [4] : 
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, (9) 



where k = y/4mE 2 /2>. From Eqs. (8) and (9), and for 
E = 0, the matrix elements of are given by: 

(^|T 4 (0)|0d> = t(0)%1) + 2r (-3qf/(Am)) h( qi ) , (10) 

where the s— wave function h{q) is the solution of 



h(q) = 



x In 
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q 2 + q' 2 — qq' ji 2 
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q k — q' + ie 

(11) 



The normalized two-body bound-state wave function, 
in the zero-range model, to be introduced in Eq.(7), is 
given by 



tt V a (h/a) 2 + p 2 



(12) 



By considering the above equations, we obtain the final 
form of the recombination parameter: 
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The numerical results for the recombination parameter 
are obtained from the solution of Eq. (11), for different 
values of fi. When \i — > oo, the results approach the scal- 
ing limit [18,19]. Therefore, the theoretical results for a 
are shown in Fig. 1 as a function of the ratio \jE 2 jE 3 . 
The calculations were performed in dimensionless units, 
such that all the momentum variables were rescaled in 
units of /i (in other words, fi = 1 in our calculations). 
So, the two-atom binding energy is decreased in respect 
to this scale. In that sense, the Thomas-Efimov states 
appear for E 2 /E li going towards zero, which is equiva- 
lent of having E 2 fixed and /j — > oo. The parameter a, 
shown in Figure 1, is obtained as a function of the most 
excited trimer state. We have performed numerical cal- 
culations with at most three Efimov states. The full cir- 
cles show the results when exists only one bound state. 
When E 2 /E^ allows two Efimov states, the results are 
represented by the solid curve, which is plotted against 
the energy of the excited state. With full squares we rep- 
resent the results when E 2 /E ll allows three Efimov states. 



The scaling limit is well approached in our calculations. 
The maximum a occurs at the threshold (E 3 = E 2 ) and 
when (E 3 /E 2 )i =0.38 [19]. So, according to this figure 
one obtains that 1 < m(a/Ti) 2 E^ < 6.9, a range consis- 
tent with Rcfs. [18,19,21]. The scaling limit for a has 
been obtained in Refs. [7,9], but without reference to the 
weakly bound triatomic molecular state. 
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FIG. 1. The dimensionless recombination parameter a as 
a function of the ratio between the binding energies of the 
diatomic and triatomic molecules. Theoretical results: full 
circles (one triatomic bound state), solid line (two triatomic 
bound states) and full squares (three triatomic states). The 
lines indicate the center of experimental data, given in Table 
I, for 23 Na, 87 Rb and 85 Rb. In case of S5 Rb we subtracted 
the contribution of the deep bound state, that was reported 
in [14]. 

In Fig. 1, we represent with horizontal lines the center 
of experimental values of a, as given in Table I, for the hy- 
perfine states 23 Na|l,-l), 87 Rb|l,-l) and 85 Rb|2,-2). 
• Using the measured values of a exp t one obtains, from 
the universal scaling plot, two weakly bound triatomic 
molecular states, denoted by E3 and E' 3 , consistent with 
these values. Considering the center of experimental val- 
ues of a, our predicted values for E 3 and E' 3 are given in 
Table I in mK. We are also giving in the table the cor- 
responding known values of the scattering lengths. The 
range for the predicted values can be easily estimated 
from Fig. 1, considering the corresponding error bars in 
a exp t- See also Ref. [22], for a recent experiment with 
ultracold thermal gas of 133 Cs|3,3), where the obtained 
values of a max , considering their systematic error limits, 
are in good agreement with theory. 

One observe that the trap dilutcness parameter is 
smaller than one in all the cases. For 85 Rb|2,— 2), we 
study a case corresponding to K 3 w 3.5 ± 1.5 x 10~ 23 
cm 6 /s, extracted from figure 2c of Ref. [13], measured 
in an ultracold non-condensed gas with external field 
B =156 G. This value of B corresponds to a = 4000a 
(ao is the Bohr radius) (See Claussen et al. [1]). As the 
resulting value of a is quite small for 85 Rb|2,— 2), one 
should expect a more significant contribution from the 
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deep bound state. Thus, we found instructive to subtract 
such contribution from a exp , which is about one unit, as 
found in Ref. [14]. However, the resulting effect in the 
determination of the trimer energy is not so dramatic, as 
seen in Fig. 1. The experimental value of a for 87 Rb|2, 2) 



does not appear in the figure, as it is well above the maxi- 
mum. By increasing the value of a from 5.8 nm to 6.8 nm 
we can make the experimental value consistent with our 
scaling limit approach. We also point out that the trimer 
can only support E 3 or E' z , not both simultaneously [19]. 



TABLE I. For the atomic species A Z\F, tuf), given in the 1st column, we present in the 6th and 7th columns our predicted 
trimer binding energies, in respect to the threshold, S3 = (E3 — E2) and S3 = (E' 3 — E2), considering the central values of 
the experimental dimensionless recombination parameters a exp t (given in the 4th column). It is also shown the corresponding 
two-body scattering lengths a (2nd column), the diluteness parameters pa 3 (3rd column), and the dimer binding energies E2 
(5th column). For 87 Rb|l,— 1), the recombination process was obtained in Ref. [11] for noncondensed (*) and condensed (*) 



trapped atoms. 


A Z\F,m F ) 


a(nm) 


pa 3 


Q-expt 


E 2 (mK) 


S 3 (mK) 


Ss(mK) 


^Na|l,-1) 


2.75 


6xl0" b 


42±12 [10] 


2.85 


4.9 


0.21 


87 Rb|l,-l> 


5.8 


lxKT 5 


52±22* [11] 


0.17 


0.39 


0.005 


87 Rb|l,-l> 


5.8 


lxKT 4 


41±17 f [11] 


0.17 


0.30 


0.013 


87 Rb|2,2) 


5.8 


4xl(T 5 


130T36 [12] 


0.17 






85 Rb|2,-2) 


211.6 


0.5 


7.84T3.4 [1,13] 


1.3xl0~ 4 


1.14xl0~ 4 


3.8 xl0~ 5 



In our predictions for the trimer's energies, except for 
85 Rb, we have disregarded the possible much smaller 
contribution of the recombination rate into deep bound 
states for a > 0, considering only recombinations into 
shallow states. When the recombination into deep bound 
states is taken into account, the curve in Fig. 1 is moved 
upward by an unknown amount. But, using the value 
found in Ref. [14], this contribution hardly is going to 
affect the extracted values for the trimer's binding ener- 
gies, given in Table I. It seems natural that, if one were 
to measure the recombination rate as a function of an 
applied magnetic field, leading to a Feshbach resonance, 
one perhaps could be able to fix this additional contribu- 
tion and determine the trimer binding energy [23]. This 
additional contribution may help to explain part of the 
measured value of a for 87 Rb|2, 2). 

In summary, in the present work, we derived the scal- 
ing dependence of the recombination parameter as a func- 
tion of the ratio between the energies of the atomic dimer 
and the most excited trimer states. The scaling func- 
tion tends to a universal function in the limit of zero- 
range interaction or infinite scattering length. The max- 
imum of the recombination rate comes at the threshold 
for the appearance of a bound triatomic molecule. In 
the cases of diluted gases of 23 Na|l, -1), 87 Rb|l, -1) and 
85 Rb|2, —2), we use the scaling function, with the corre- 
sponding known experimental values of the recombina- 
tion rates and two-atom scattering lengths, to predict 
for the first time the binding energies of weakly bound 
trimers in ultracold traps. We stress that the possi- 
ble contribution of a deep bound state in our predic- 
tions is expected to be relatively small, as verified for 
85 Rb|2, —2), in a particular trap. We also note that for 
the 85 Rb|2,— 2) the diluteness parameter is about 0.5, a 



value that may be considered near the limit of validity 
of the present approach, which does not include higher 
order correlations between the particles. 

Finally, we would like to remark that, at a first sight, 
one could think that formation of trimers requires four- 
body collisions, which are very unlikely unless the density 
is high. However, the recent experimental results given 
in Ref. [1] indicate formation of molecules in the trap, 
as also discussed in Ref. [24]. Therefore, other collision 
processes like dimer-dimer or dimer and two atoms could 
also lead to trimer formations, enhancing the possibility 
of producing trimers in a trapped ultracold gas. 
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